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ABSTRACT 

We explore the ability of weak lensing surveys to locate massive clusters. We use both analytic 
models of dark matter halos and mock weak lensing surveys generated from a large cosmo- 
logical TV-body simulation. The analytic models describe average properties of weak lensing 
halos and predict the number counts, enabling us to compute an effective survey selection 
function. We argue that the detectability of massive halos depends not only on the halo mass 
but also strongly on redshift at which the halo is located. We test the model prediction for 
the peak number counts in weak lensing mass maps against the mock numerical data, and find 
that the noise due to intrinsic galaxy ellipticities causes a systematic effect which increases the 
peak counts. We develop a correction scheme for the systematic effect in an empirical man- 
ner, and show that, after the correction, the model prediction agrees well with the mock data. 
The mock data is also used to examine the completeness and efficiency of the weak lensing 
halo search with fully taking into account the noise and the projection effect by large-scale 
structures. We show that the detection threshold of S/N = 4^5 gives an optimal balance 
between completeness and efficiency. Our results suggest that, for a weak lensing survey with 
a galaxy number density of n g = 30 arcmin~ 2 with a mean redshift z = 1, the mean number 
of halos which are expected to cause lensing signals above S/N = 4 is iVh a io(S/N > 4) = 37 
per 10 deg 2 , whereas 23 of the halos are actually detected with S/N > 4, giving the effective 
completeness as good as 63%. On the other hand, the mean number of peaks in the same area 
is iVp ea k = 62 for a detection threshold S /N = 4. Among the 62 peaks, 23 are due to halos 
with the expected peak height S/N > 4, 13 are due to halos with 3 < S/N < 4 and the 
remaining 26 peaks are either the false peaks due to the noise or halos with a lower expected 
peak height. Therefore the contamination rate is 42% (this could be an overestimation). Weak 
lensing surveys thus provide a reasonably efficient way to search for massive clusters. 

Key words: cosmology: theory — gravitational lensing — galaxies: halos — dark matter — 
large-scale structure of universe 



1 INTRODUCTION 

Recently, it has become feasible to locate massive clusters di- 
rectly as density enhancements using weak gravitational lensing. 
Miyazaki et al. (2002) indeed discovered many clusters in a weak 
lensing halo survey over a 2. 1 square degree field, proving the abil- 
ity of weak lensing to identify massive clusters (see also Wittman 
et al. 2001). Unlike conventional optical or X-ray selected cluster 
catalogs, weak lensing cluster catalogs are free from a bias toward 
the luminous objects because the cluster finding is not based on the 
flux enhancement but on the projected mass enhancement. This is 
a great advantage of gravitational lensing. However, lensing has its 
own disadvantages as every other cluster survey technique (opti- 
cal, X-ray and Sunyaev-Zel'dovich survey) does; for example, the 
projection effect from unrelated structures in the same line-of-sight 
(Reblinsky & Bartelmann 1999; Metzler, White & Loken 2001). 



Since clusters of galaxies are one of most important cosmologi- 
cal probes, it is of fundamental importance to construct large unbi- 
ased samples. Given advantages and disadvantages in each survey 
technique, it is desirable to combine, in a complementary manner, 
several techniques to construct an unbiased catalog. Such catalogs 
provide a valuable data set to investigate, for example, the physical 
state of the intra-cluster medium. To do this, it is important to cor- 
rectly understand the selection function and completeness of each 
survey technique. 

It has been often argued that weak lensing provides a truly 
mass-selected sample of clusters. This is not strictly true; weak 
lensing measures the two-dimensional tidal field which is a line- 
of-sight projection of the three-dimensional tidal field weighted by 
distance ratio DiDi s /D s (see §2 for definitions and details; see 
Mellier 1999; Bartelmann & Schneider 2001 for reviews of weak 
lensing). In simpler terms, the amplitude of the lensing signal from 
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a cluster is not solely determined by its mass, but depends also on 
the redshift and the shape of the gravitational potential due to the 
cluster. In addition, both foreground and background large-scale 
structures contribute to the lensing signal, and this projection ef- 
fect usually adds a noise to the lensing signal from a halo. Further, 
since the weak lensing signal is measured from tiny coherent dis- 
tortions in galaxy shapes, intrinsic ellipticities of galaxies introduce 
an irremovable noise. Thus, the detectability of weak lensing halo 
depends on all these factors. 

The primary purpose of the present paper is to explore the 
ability of weak lensing surveys to locate massive clusters. We es- 
pecially address the following three points; (1) to examine the se- 
lection function of weak lensing cluster surveys, (2) to develop a 
theoretical model that describes the weak lensing halo counts, and 
(3) to examine how the detailed structure of clusters, projection ef- 
fect and noise affect weak lensing cluster surveys under a typical 
observational condition from a ground based 4 — 10m class tele- 
scope. The last point is especially important to understand a bias 
that weak lensing halo catalogs may have. To address these things, 
we use both simple analytic descriptions of dark matter halos and 
mock numerical weak lensing survey data. The former offers a use- 
ful way to compute expected lensing properties of massive halos. 
We use them to compute the selection function. We also develop a 
theoretical model for the number counts of the weak lensing halos 
based on the analytic descriptions. On the other hand, the latter al- 
lows us to examine all the factors listed above in a direct manner. 
Mock numerical catalogs are generated using weak lensing ray- 
tracing simulations. We extensively use the mock data and the halo 
catalogs directly produced from the simulation outputs to examine 
completeness and efficiency of the weak lensing halo survey. The 
mock data is also used to test the model prediction of the weak 
lensing halo counts. 

White, Van Waerbeke, & Mackey (2002) and Padmanabhan, 
Seljak & Pen (2003) addressed the ability of weak lensing surveys 
focusing on the projection effect. In those studies, the complete- 
ness of weak lensing surveys is defined as the fraction of detected 
halos with mass above a certain value relative to halos that lie in 
a given volume. This definition compromises two completely dif- 
ferent elements, (i) the selection effect and (ii) the effects from in- 
dividuality of the halo mass distribution, the projection and noise. 
For our purpose, these elements needs to be examined separately. 
We first compute the selection function using the analytical descrip- 
tions of dark matter halos which enables us to define "potentially 
detectable halos". Here, "potentially" refers to an ideal case in the 
absence of noise. Then, we shall adopt the different definition of the 
completeness that is the fraction of actually detected weak lensing 
halos relative to the potentially detectable halos. 

Our results have important implications to real observations. 
When making an observational strategy, one would like to com- 
pute the expected number of detectable clusters using, for exam- 
ple, a simplified model, and also needs to estimate the expected 
completeness and false detection rate. Our theoretical model of the 
weak lensing halo counts may be useful for the former case, and 
our results from the mock data are directly applicable to derive 
the latter estimates. Also, from the observational point of view, we 
adopt a directly observable quantity, a peak height of weak lensing 
mass map divided by the noise root-mean-square (RMS), that is 
the signal-to-noise ratio, S/N, of the detection, as the fundamental 
estimate of the lensing signal. 

The rest of this paper is organized as follows. In §2, we sum- 
marize lensing properties of massive halos assuming the univer- 
sal density profile proposed by Navarro, Frenk & White (Navarro, 



Frenk & White 1996; 1997). Using the expected properties, we 
compute the selection function of the weak lensing cluster survey. 
Then, we combine the detectability of weak lensing halos with the 
halo mass function based on the Press-Schechter approach (Press 
& Schechter 1974; we adopt the modified fitting function by Sheth 
& Tormen 1999) to compute the weak lensing halo counts and 
the selection functions with respect to the halo mass and redshift. 
This simple analytical approach should offer a theoretical basis to 
explore properties of weak lensing halos. In §3, we describe our 
numerical experiment of weak lensing surveys and weak lensing 
ray-tracing technique. Statistical properties of halo catalogs are ob- 
tained and compared with the theoretical predictions derived in §2. 
Using the mock weak lensing survey data, in §4, we examine the 
number counts of peaks in the weak lensing convergence map, and 
the theoretical prediction developed in §2 is tested. In §5, we in- 
vestigate correspondences between halos in the halo catalog and 
peaks identified in the lensing convergence map. In particular, we 
examine the completeness and efficiency of the weak lensing halo 
search paying special attention to the points (1-3) above. Summary 
and discussion are given in §6. In Appendix A, we present some 
lensing convergence maps in which a missing halo or a false peak 
exists for illustrative examples of irregular systems. In Appendix 
B, we develop a correction scheme to the theoretical prediction of 
the weak lensing peak counts. In Appendix C, relations between 
the halo shape (and orientation) and the strength of lensing signal 
are examined. 

Throughout this paper, we work with the standard flat ACDM 
cosmological model with the density parameter fl m — 0.3, the 
cosmological constant £7a = 0.7, and Hubble constant Ho = 
lOO/i kms _1 Mpc _1 with h = 0.7. We adopt the fitting func- 
tion of the CDM power spectrum by Bardeen et al. (1986) with the 
normalization of ag = 0.9. We take the observational parameters 
of a typical weak lensing survey from a ground based 4 — 10m 
telescope at the present. To be specific, we consider an imaging ob- 
servation with a limiting magnitude fainter than R — 25.5 mag in a 
sub-arcsecond seeing condition, which provides the galaxy number 
density n g > 30 arcmin - with the mean redshift z ~ 1. 



2 PREDICTIONS FROM STANDARD MODELS 
2.1 Lensing properties of an NFW halo 

The universal density profile of dark matter halos proposed by 
Navarro, Frenk & White (1996; 1997; NFW hereafter) is given by 



PnfwW 



x(l + x) 2 ' 



(1) 



where r s is the scale radius. It is convenient to introduce the con- 
centration parameter c n f w = r v i r /r s where r v i r is the virial radius. 
Bullock et al. (2001) found from iV-body simulations that the con- 
centration parameter is related to halo mass as 



Cnfw(M, z) = 



M 



1 + z V IQ^h-^Mp 



(2) 



where we set c* = 8 for the ACDM model. It is important to notice 
that there is a relatively large scatter in this relation (Jing 2000). The 
mass enclosed within a sphere of radius r v i r gives the virial mass 
by definition: 



Mvir 



l0g(l + C nfw ) 



Cnfw 



1 + Cnf 



(3) 



© 2003 RAS, MNRAS Q0Q.ITTI221 



Searching for massive clusters in weak lensing surveys 3 



The virial mass is also defined by the spherical top-hat collapse 
model as 



Mvir = -^-<5vir(2)P0JVi, 



(4) 



where 5 v i r is the threshold over-density for spherical collapse (see 
Nakamura & Suto 1997 and Henry 2000 for useful fitting func- 
tions). From equations J3} and J4} one can express p s in terms of 
8 v ii(z) and c n f w . Introducing S s = p s /po — 1, one finds, 



5s = 



r 3 

u nfw 



3 log(l + 

^nfw) Cnfw/ ( 1 "T" C n f w ) 



(5) 



Although the NFW profile formally extends to the infinity, an 
actual halo has a finite extent. We consider a truncated NFW halo, 
the density profile of which is truncated at the virial radius. The sur- 
face mass density of the truncated NFW halo is given by the line- 
of-sight integration of the density profile (Takada & Jain 2003): 



E(y) 

with 



dz p ntw (y,z) = 2p s r s f(y), y=— ,(6) 



" (l-W^Xl + Cnfw) + (1-1/2)3/2 



r arccosh 



V + e nfw 
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for y < 1 



/(v) = < 



V^fw-" 2 _ 1 
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(7) 



(l-y2)(l + c nfw ) (!/ 2_ 1) 3/2"— j/(l + c nfw ) ' 

for 1 < y < c, 

0, for y > c. 

Note that the above expression differs from that by Bartelmann 
(1996) and Wright & Brainerd (2000) which is for the non- 
truncated NFW halo and gives a larger surface mass density than 
ours. The lensing convergence profile is given by 



<V) = ^ = «./(»), 

with the critical surface mass density 

" 2 D s 



4itG DiDi. 



and 
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where a is the scale factor normalized at the present time, D de- 
notes the angular diameter distance, x denotes the comoving radial 
distance (which is equal to the comoving angular diameter distance 
for a flat universe), and the subscripts I and s stand for lens and 
source, respectively. It is worth noting that k scales with c n f w ap- 
proximately as k oc c n f w for 1 < Cnfw 5 and k oc c nfw for 
5 < Cnfw i$ 30. This clearly proves that the convergence ampli- 
tude is sensitive to the halo concentration. 

For the number density of source galaxies and the intrinsic 
ellipticity RMS we are interested in, a raw (un-smoothed) conver- 
gence map is very noisy. In order to have a lensing mass map with 
reasonably small amplitude of the noise, we smooth the conver- 
gence map with a Gaussian window function, 



W G {x) 



1 



(ii) 



CD 

o 




Figure 1. The solid line shows the lensing signal from an NFW halo 
smoothed with the Gaussian kernel (eq. 1131 ) against the smoothing length. 
For the concentration parameter, c ntw = 5 is adopted here. The dotted line 
shows scaling of the noise RMS with the smoothing scale (<r no ; se oc 9^, , 
see eq. 1121 ). 



The Gaussian smoothing is indeed adopted in actual weak lensing 
surveys (e.g., Miyazaki et al. 2002). As shown in Van Waerbeke 
(2000), if the intrinsic ellipticities are uncorrelated between differ- 
ent source galaxies, statistical properties of the resulting noise field 
can be described by the Gaussian random field theory (Bardeen et 
al. 1986; Bond & Efstathiou 1987), on scales where the discrete- 
ness effect of source galaxies can be ignored. The Gaussian field is 
specified by the variance of the noise field that is specified by the 
number of galaxies contained within a smoothing aperture (Kaiser 
& Squires 1993; Van Waerbeke 2000): 

-2 -i 



2 2ir0^n 9 ! 



(12) 



where <r £ is the RMS amplitude of the intrinsic ellipticity distri- 
bution and n 9 is the number density of source galaxies. Through- 
out this paper we adopt a e — 0.4 and n g = 30 arcmin~ 2 , which 
are typical values for a survey with a limiting magnitude fainter 
than R = 25.5 mag in a sub-arcsecond seeing condition. With this 
choice, the noise RMS (J no i sc = 0.02 x (1 acrmin/^c). It is worth 
noting that er no i sc is a direct observable in sense that it can be esti- 
mated from a reconstructed weak lensing mass map obtained after 
randomizing the orientation of source galaxy images. 

We need to also take into account the smoothing effect on the 
lensing signal from a halo. The peak amplitude of the convergence 
profile of halos is modified by smoothing into 



kg(9g) = / d 2 <pW G (<f>;e G )K(<f>), 



(13) 



where the center of the smoothing kernel is set to the halo center. 
Hereafter we refer to this peak amplitude as "the lensing signal 
due to a halo". Figure Q plots the lensing signal as a function of 
the smoothing scale. The dotted line shows a scaling of the noise 
RMS with the smoothing scale (a no i sc oc 6q ). The comparison 
between the solid and dotted curves thus shows how the signal- 
to-noise ratio (S/N) varies with the smoothing scale. Clearly, an 
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Figure 2. The gray scale with contour lines show the angular scale radius 
B in unit of arcminutes, in the halo mass (Af^aio) and redshift (zr) plane. 
The source redshift is fixed to be z s = 1. 




0.5 1 



Figure 3. The gray scale with contour lines shows the S/N value for weak 
lensing halo detection defined by eq. 1141 . The S/N value is computed 
assuming the universal NFW profile. We adopt, throughout this paper, 
u t = 0.4, n g = 30arcmin -1 and 9q = 1 arcmin. 

optimal S/N can be attained for 9c — (1 — 2) x 8 S . Figure[2]shows 
the contour plot of 8 S in the halo mass and redshift plane, where 
the mean mass-concentration relation given by eq. J2j is employed. 
Note that throughout this paper, the source redshift is fixed to be 
z s — 1, which is a typical value for the mean redshift of source 
galaxies for an actual survey with limiting magnitude R ~ 25.5. 
Combining Figures □ andU it is suggested that dc ~ 1 arcmin 
would be an optimal choice for efficient survey of massive halos 
with Mhato > 10 14 /i _1 Mq at intermediate redshifts. 

We define the signal-to-noise ratio, S/N, for a weak lensing 
halo detection by 1 



1 Note on the notation issue: Throughout this paper, we denote the S/N of 
weak lensing signals by v and distinguish its type by the subscript; f n f w 
stands for the S/N computed assuming the universal NFW profile, while 
Vl ens and f no i S y for the S/N measured in the noise-free and noisy realiza- 
tions of numerical weak lensing mock data, respectively (see J3]and JIJ. 




0.5 



Figure 4. The gray scale with solid contour lines shows the num- 
ber density of halos per 1 deg 2 area, a given redshift interval dz and 
a given logarithmic mass interval dlogM: d 2 N^^/ dz / d\og M = 
dV I dz dn^ic/d log M . This is the integrand of eq. 1151 . The dashed con- 
tours show the expected peak height u n{v/ from 3 to 9 at an interval of 
A^nfw = 1. Integrating this number over the peak height above a thresh- 
old value v t h gives the expected number of peaks per 1 deg 2 as given by 
eq. (HI. 



V noise 

Figure|3|shows the contour plot of the S/N value computed assum- 
ing the universal NFW profile in the halo mass and redshift plane. It 
is evidently seen from this figure that the lensing signal not only de- 
pends on the mass but also on the redshift. This redshift dependence 
enters mainly through the lensing efficiency function DiDi 3 /D s 
(see eqs. (9) and 1101 ). This makes the mass sensitivity of weak 
lensing surveys strongly redshift dependent. In other words, the 
lower mass limit of the detectable halos depends strongly on the 
redshift. This figure also implies that the effective survey volume of 
a weak lensing halo search depends on the minimum halo mass that 
one aims at locating. 

It should be noticed that we adopt the mean mass- 
concentration relation to compute the model prediction for the S/N 
value, though the relation has large scatters (Bullock et al. 2001; 
Jing 2000). Since the lensing signal, kg, depends on the concen- 
tration parameter as kg oc cz^ , the scatter likely translates to 
scatter in the kg value or equivalently in the S/N. It should be also 
noted that the lensing signal for a realistic halo is affected by sub- 
structures and asphericity of the mass distribution seen in iV-body 
simulations. We will carefully investigate these effects, using mock 
data of weak lensing survey. 

2.2 Model prediction of weak lensing halo counts 

Halos that induce a sufficiently high lensing signal are likely to 
be identified as peaks in the convergence map, even if the map is 
contaminated by the noise due to intrinsic ellipticities. Hence, it 
is convenient to work with the number counts of peaks above a 
certain threshold value: N(u > !/ tn ; #g)- In order to compute the 
number counts, we assume that a high peak is generated by the 

We will also use the notations u t ^ and u mln to denote a threshold value of 
the detection and the minimum value to be set (see [|l]and J5J, respectively. 
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Table 1. Summary of the weak lensing peak counts computed using eq. H5\ 
for various smoothing scales. Observational parameters adopted are cr e = 
0.4, n g = 30arcmin -1 , and a fixed source redshift of z s = 1. The second 
and third columns are for the counts v > 5 and v > 4, respectively. 



9q [arcmin] N(u > 5) [lsq. 


deg" 1 ] N{u > 4) [lsq. deg- 1 ] 


0.5 


1.47 


3.38 


1 


2.08 


4.05 


2 


1.76 


3.15 


3 


1.30 


2.27 


4 


0.97 


1.67 


Table 2. Summary of the 


angular number density (per 1 deg 2 ) of halos with 


mass more massive than 


a given 


minimum halo mass (the first column). 


The second and third columns are 


for the maximum redshifts z < 1 and 


2 < 0.7, respectively. 






Halo mass [/i _1 JW Q ] 




[lsq. deg -1 ] 




2 < 1 


2 < 0.7 


M > 5 x 10 13 


55 


30 


M > 1 x 10 14 


17 


10 


M > 2 x 10 14 


4.1 


2.7 


M > 5 x 10 14 


0.36 


0.28 



lensing signal from a single massive halo and that the number of 
false peaks due to noise (i.e., contamination) is much smaller than 
that of real peaks, as will be shown below in detail. For our fiducial 
choice of the smoothing scale 9q = 1 arcmin, this assumption is 
reasonably valid for sufficiently high peaks (say v > 4), and thus 
we focus on such high peaks. In this case, the number density of 
peaks with height greater than a threshold v x \^ per unit solid angle is 
given by (Kruse & Schneider 2000; Bartelmann, King & Schneider 
2001) 



N (v > v t h-, 0G,z s ) = 



dx 



dV 



dM 



drihaio(M, z{x)) 



dx J dM 
xH(y{M,z) - i/ t h), 



(15) 



where \ ls the radial comoving distance, dV/dx is the comoving 
volume element per unit solid angle (dV/dx = X 2 f° r a A at cos ~ 
mology), drihaio/dM is the halo mass function, for which we adopt 
the fitting function by Sheth & Tormen (1999), and TL(x) is the 
Heaviside step function (TL(x) = 1 for x > and otherwise). 
The differential counts are obtained by 



n p eak(^) = 



dN(u > u th ; 9 G , z s ) 



dv 



(16) 



Table □ summarizes the predicted weak lensing peak counts 
for various smoothing scales. It is shown that a choice of Qq — 1 ar- 
cmin gives the largest counts, for cr e = 0.4 and n g = 30 arcmin -2 . 
In the following discussion, we will adopt this smoothing scale. 

In Figure [4] we plot the selection function of weak lensing 
halo search with respect to the halo mass and redshift. The gray 
scale with solid line contours shows the number density of halos 
per 1 deg 2 against a given redshift interval of dz and a given loga- 
rithmic mass interval of dlog M: 



d 2 N h 



dV dWhalo 



(17) 

dz dlogM dz dlogM 

Note that this is the integrand of eq. 1151 . It is clear that, for a given 
halo mass, the number density rapidly increases with redshift at 



10 2 



10 



-a 1 

T3 



0.1 



M (l =5xl0 13 : 



lxlO 14 




Figure 5. The solid lines show the redshift distribution of halos with a peak 
height above a threshold value u t h per 1 deg 2 , dN(u n [ v > u t ^)/dz/l 
deg 2 , computed using eq. 1151 . The upper and lower solid lines are for 
= 4 and 5, respectively. Note that we employed 9q = 1 arcmin and 
z s = 1. The broken lines show the same redshift distributions but for ha- 
los with masses above a given minimum mass, dN^i^M > M t ^)/dz/l 
deg 2 . 



low redshifts z < 0.2, and then gradually decreases toward higher 
redshifts. This is because the volume element, dV oc x 2 ^X' i n " 
creases with redshift, whereas the halo number density decreases. 
The dashed contours show the expected peak height ^ n f w for the 
universal NFW halo at intervals of A!/ n f w = 1. The comparison 
between the solid and dashed contours clarifies that, if threshold 
peak height i^ t h = 4 — 5 is employed, the weak lensing halo search 
is most sensitive to halos that are located at the redshift between 
0-1 ^ z < 0.7 and have masses greater than 10 14 h" 1 Mq. 

Let us look further into the redshift- and mass-dependence of 
the weak lensing selected halo counts. The solid curves in Figure 
|5|show the redshift distribution of halos selected by the weak lens- 
ing signal. These curves can be compared with the distribution of 
halos above a certain mass threshold (broken lines). A remarkable 
difference between these two is evidently seen at redshifts z > 0.5; 
the distribution of weak lensing selected halos drops very rapidly, 
while the mass selected ones decrease very gradually. This is again 
because the lensing efficiency peaks around z = 0.4 for source 
galaxies of z = 1. Figure|6|shows essentially the same information 
but from a different angle. The left panel shows the mass function 
of halos with peak height exceeding a threshold value (the cases 
of ^ t h = 4 and 5 are displayed). The lower mass limit in the mass 
function is determined by the detection threshold, while the decline 
in the high mass side (which has an approximate exponential slope) 
is due to the property of the halo mass function (i.e., the exponen- 
tial shape of its high mass end, see Kruse & Schneider 2000). Right 
panels of Figure|6|shows the mass distribution of halos for different 
redshift intervals as labeled in each plot. The solid (dashed) lines 
are for halos with u n f v > 5 (z^ n fw > 4), while the dotted lines 
are for all the halos in the corresponding volume. The upper two 
plots clearly show that, up to z = 0.4, almost all the halos with 
M > 2 x lO 14 ft _1 M0 indeed produce high peaks with u n f w > 5. 
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10 14 10 15 10 14 10 15 



M hal0 [h'-M e ] M halo [h'-M ] 

Figure 6. Left Panel: Mass function of halos with u n f w > u t ^ computed 
using eq. 1151 - The threshold values are u t h = 4 and 5 for the dashed 
and solid lines, respectively, The source redshift is taken to be 2 S = 1. 
Right Panels: Same as the left panel but divided into four ranges of the halo 
redshift as denoted in each plot. The dashed lines show the mass function of 
halos within the corresponding volume (which is equivalent to set z/ t h > 0, 
that is, no limit on the weak lensing signal is imposed). Note that the scales 
of y-axis are labeled on the right side of the plots. 



However, for higher redshifts, only very massive halos can produce 
high peaks. 

To summarize the results shown in Figures|4]-|6| we argue that 
weak lensing is not an efficient technique of finding high-redshift 
halos (say z > 0.6), except for very massive halos of M > 5 x 
lO 14 /i _1 M ,if a mean source redshift is around z 3 = 1. It should 
be noticed that the ability of the weak lensing halo surveys depends 
on observational conditions. As a survey is deeper, the halo search 
more enables us to detect less massive halos and halos at higher 
redshift. 

Table [2] gives the number density of halos with mass above 
a given minimum mass, which is obtained by integrating the halo 
mass function up to a given maximum redshift. Comparison with 
Table Q provides a rough estimate of how weak lensing survey 
misses mass-selected halos in the given survey volume. A large dis- 
crepancy between the halo counts and the weak lensing halo counts 
is seen, which arises from the redshift dependence of the lensing 
selection function as shown in Figures|5]and|6| 



3 NUMERICAL SIMULATION AND MOCK CATALOGS 

The theoretical model prediction discussed in the previous section 
provides us with a basic selection function of the weak lensing clus- 
ter survey. However, since our model is based on the simple aver- 
aged descriptions of the dark matter halos, it does not take into 
account the individuality of halos which causes a substantial scat- 
ter in the lensing signal and may even cause a bias. Furthermore, 
there are two additional sources of scatter; the projection effect and 
the intrinsic ellipticity noise. In this section, we will examine these 
effects in detail using mock numerical experiments. 



3.1 TV-body simulation 

We used TV-body data from Very Large TV-body Simulation (VLS) 
carried out by the Virgo Consortium (Jenkins et al. 2001, and see 
also Yoshida et al. 2001 for simulation details). The simulation 
was performed using a parallel P 3 M code (MacFarland et al. 1998) 
with a force softening length of l EO { t ~ 30 h~ kpc. The simula- 
tion employed 512 3 CDM particles in a cubic box of 479 /i _1 Mpc 
on a side. It uses a flat cosmological model with Qo = 0.3, 
Qa = 0.7, and h = 0.7. The initial matter power spectrum was 
computed using CMBFAST (Seljak & Zaldarriaga 1996) assum- 
ing a baryonic matter density of Q,h = 0.04. The normalization 
of the power spectrum is set by as = 0.9. The particle mass 
(wipart = 6.86 x 10 10 /i _1 Mq) of the simulation is sufficiently 
small to guarantee practically no discreteness effect on dark-matter 
clustering on scales down to the softening length in the redshift 
range of interest for our purposes (Hamana, Yoshida & Suto 2002). 



3.2 Weak lensing ray-tracing simulation 

The multiple-lens plane ray-tracing algorithm we used is detailed in 
Hamana & Mellier (2001; see also Bartelmann & Schneider 1992, 
Jain, Seljak & White 2000 and Vale & White 2003 for the the- 
oretical basics and technical issues); here we thus describe only 
aspects specific to the VLS TV-body data. In order to generate the 
density field between 2 = and z ~ 1, we use a stack of 6 snap- 
shot outputs from two runs of the TV-body simulation, which differ 
only in the realization of the initial fluctuation field. Each cubic 
box is divided into 4 sub-boxes of 479 2 x 119.75/i~ 3 Mpc 3 with 
the shorter box side being aligned with the line-of-sight direction. 
TV-body particles in each sub-box are projected onto the plane per- 
pendicular to the shorter box side and thus to the line-of-sight di- 
rection. In this way, the particle distribution between the observer 
and 2 ~ 1 is projected onto 19 lens planes separated equally by 
119.75 ft _1 Mpc. Note that in order to minimize the difference in 
redshift between a lens plane and an output of TV-body data, only 
one half of the outputs (i.e. two sub-boxes) from 2 = TV-body 
box was used. 

The particle distribution on each plane is converted into the 
surface density field on a 4096 2 regular grid using the triangu- 
lar shaped cloud (TSC) assignment scheme (Hockney & Eastwood 
1988). 4096 2 grid is chosen to maintain the resolution provided 
by the TV-body simulation as well as to remove at the same time 
the shot noise due to discreteness in the TV-body simulation. The 
computation follows the procedure described in Hamana & Mellier 
(2001) and Jain et al. (2000). 

Having produced surface density fields on all lens planes, we 
trace 512 2 rays backwards from the observer's point to source plane 
using the multiple-lens plane algorithm (e.g. Schneider, Ehlers & 
Falco 1992). The initial ray directions are set on 512 2 grids with 
a grid spacing of 0.24 arcmin, and thus the total area covered by 
the rays is 2.048 2 square degrees. We made 50 realizations of the 
underlying density field by randomly shifting the simulation boxes 
in the direction perpendicular to the line-of-sight using the periodic 
boundary conditions of the simulation 2 . Note that the lens planes 
coming from the same box are shifted in the same way in order to 
maintain the clustering of matter in one box. 

In the following analyses, we use the lensing convergence map 

2 The weak lensing ray-tracing simulation data used in this paper are avail- 
able from T. Hamana on request, hamanatk@cc.nao.ac.jp. 
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Figure 7. The number counts of the lensing peaks above a given threshold 
^nfw P er a 4 square degrees area. The solid curve shows the model predic- 
tion computed from eq. 1151 , while the filled circles with error bars denote 
the mean and RMS among 50 realizations of mock simulation data. Note 
that the peak height of a halo is estimated from the halo mass and redshift 
in the manner described in the text ( A3. 41 . 



of the source plane closest to z — 1. This is, we call, the noise-free 
lensing map. 



3.3 Adding noise 

To simulate a more realistic lensing survey, we add noises due to in- 
trinsic ellipticities of source galaxies to the simulated lensing map. 
Actual observations are further affected by other systematics such 
as imperfect conditions of telescope and/or of electric devices of 
detector. Examining these effects is beyond the scope of the present 
paper, and thus we do not consider such observational effects. 

As extensively shown in Van Waerbeke (2000), the noise field 
due to the intrinsic ellipticities is well approximated by the Gaus- 
sian random field on angular scales of our interest, if the elliptic- 
ities are uncorrelated between different sources (see also Jain & 
Van Waerbeke 2000). The Gaussian noise added to our k map has 
variance 



2 n a e: 



2 ■ 

9 w pix 



(18) 



where 



= 0.24 arcmin is the pixel size of our k maps. The 



noise RMS on pixels is a pix = 0.22 which is much larger than the 
RMS of the lensing convergence field. 



3.4 Halo catalogs on the light-cone 

We identify dark matter halos in 5 low redshift outputs of the 
iV-body simulation using the standard friends-of-friends algorithm 
with a linking parameter of b = 0.164 (in units of the mean particle 
separation). Halos with mass greater than 1 x 10 13 /i -1 M©, which 
corresponds to the mass of 146 simulation particles, are used in the 
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Figure 8. The frequency distribution of numbers of halos with the expected 
peak height (see i|3.4l for its definition) f n f w > 5 from 50 realizations of 4 
square degrees. 



analyses below. We define the halo center to be the center-of-mass 
position of member particles. 

In order to evaluate the halo shape, we compute the inertia 
tensor of the mass distribution; 



dx 3 {xi — Xi)(Xj — xj)p(x) 



N, 



(19) 



where i,j run from 1 to 3, Xi denotes the center-of-mass position 
and N p is the number of member particles of the halo. From the 
inertia tensor we compute the eigen vectors and convert them to 
the axial ratios a/c and a/b (a ^ b ^ c). We also compute the 
trace of the inertia tensor, 



Q = -(Tii +I22 + I33) 



(20) 



This is compared with the expectation value for the NFW halo, 

Qnfw 

O J 

= Mvi r rv ir g(c n fw), (21) 



— J dr r p nfw (r) 



with 

g(Cnfw) 



c 3 f 

°niw 



3c" 



6c n fw + 6(1 + )log(l + 



6c2 fw [(l + c nfw )log(l + 



(22) 



Note that q(c n f w ) can be approximated by g(c n fw) — 0.22c~ f ^ 
for 3 < c n f w < 30. The ratio Q/Qnfw gives an estimate of the 
concentration of the mass distribution. We also define the following 
quantity as an estimate of the line-of-sight elongation of the halo 
mass distribution: 

2/33 



R = 



In + I22 ' 



(23) 



© 2003 RAS, MNRAS 000.1711221 




Figure 9. Upper Panel: The scatter plot showing the relation between 
Q/Qnfw and the axial ratio a/c for halos with M > 3 X 10 13 /i — 1 Mq 
in the z = 0.1 output. Here Qnfw is defined by eq. 1211 and is computed 
using the relation eq. J4}- Lower Panel: The PDF of Q/Qnfw °f halos 
with M > 3 X 10 13 h~ 1 MQ, the solid and dashed lines for z = 0.1 and 
z = 0.5 outputs, respectively. 

where X3 is taken by the direction parallel to the line-of-sight. Al- 
though this is rather a crude estimation, R > 1 or R < 1 generally 
implies halos elongated along or perpendicular to the line-of-sight. 

Halo catalogs on the light-cone are generated by stacking the 
simulation outputs in the same manner as in the ray-tracing experi- 
ments. We use the 50 realizations which have the same underlying 
matter distribution as used in the ray-tracing data. This combined 
data set allows us to directly relate the distribution of dark matter 
halos to the weak lensing convergence map. 

In summary, each halo in the mock catalogs has data on the 
mass, the estimates of the halo shape Q and R defined by eqs. J20> 
and <23> , respectively, the redshift computed from the radial dis- 
tance to a halo, and the angular position converted from the spa- 
tial position. In addition, for each halo, we compute the virial ra- 
dius and the concentration parameter using eqs. J4j and J2j, from 
the halo mass and redshift of TV-body data output time (which are 
z = 0, 0.1, 0.3, 0.5 and 0.8). These values allow us to evaluate the 
expected lensing signal using eq. J 1 3 1 with eq. {8j, which is correct 
only if the halo has the NFW profile. Similarly, we also compute 
Qnfw using eq. J2H . Moreover, we introduce a flag parameter 
for each halo, to quantify the following possible problem that arises 
from our use of multiple-lens plane ray-tracing algorithm. As a re- 
sult of projection of TV-body particles into 4 lens planes, member 
particles of given one halo located close to the boundary could be 
separated into different planes. Such a split halo would generate an 
artificially deformed (or even separated) lensing signal in a k map. 
We assign flag = 1 for such halos whose virial radius is larger 
than the distance between the halo center and the closest plane, 
flag = 2 for halos whose scale radius is larger than the distance to 
the closest plane, and flag = otherwise. The fraction of flagged 
halos is very small, 0.88% for flag = 1 and 0.19% for flag = 2, 
thus it does not have significant influence on our analyses below, 
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R=2I 33 /(I n +I g2 ) 

Figure 10. The PDF of the R defined by eq. for halos with M > 
3 X 10 13 /i -1 Mq at z = 0.1 (solid line) and z = 0.5 (dashed line) outputs. 
Note that the median of the distributions is ~ 0.92, and 55% of the all 
samples has R < 1. 

but we carefully check if flagged halos cause any artificial eccen- 
tric signals in the lensing convergence map. 

The expected lensing signal of each halo is converted into the 
peak height in the same way as in eq. 1 141 . which we denote as 
v n fw Total numbers of halos with ;/ n f w > 5 and !/ n f w > 4 in our 
50 realizations are 459 and 871, respectively, which are sufficiently 
large to allow the statistical analyses of weak lensing halos. 

The redshift and mass distributions of halos in the halo catalog 
agree well with the model predictions in £12.21 Figure^plots the ac- 
cumulate number counts of halos against the expected peak height, 
v a tw The solid curve is the theoretical prediction of eq. 1 1 51 and 
agrees well with the simulation result. This agreement is not sur- 
prising, because ingredients of the theoretical model (the universal 
profile and halo mass function) are calibrated by TV-body simula- 
tions. 

The counts of halos with z/ n f w > 5 in 2.048 2 square degrees is 
found to be TV (^ n f w > 5) = 9.18±3.63, where the error represents 
the RMS among 50 realizations. This is in good agreement with the 
theoretical prediction of 8.71. However, the RMS is somewhat large 
and slightly larger than the Poisson fluctuation. Figure|8|shows the 
distribution of numbers of halos with i/ n f w > 5 in each realization, 
quite a wide spread is evidently seen. This is mainly due to the very 
strong clustering of massive halos. Accordingly, given a limited 
survey area, the halo counts may significantly vary from survey to 
survey. 

The lower panel of Figure [5] shows the PDFs of Q /Qnfw 
computed from massive halos with M > 3 x 10 13 Tt -1 Mq (Qnfw 
is computed using the mean mass-concentration relation eq. (2j). 
The PDF peaks at Q/Qnfw ~ 1 as expected, but has rather 
broad distribution. No significant difference in the PDFs between 
the five different redshift boxes was found. The upper panel shows 
the scatter plot in Q/QNFW-a/c plane for halos with M > 
3 x 1O 13 /i _1 M from z = 0.1 output. It is clear that the mass 
concentration defined by Q/Qnfw correlates with the axial ratio, 
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Figure 11. Upper-left: The celestial distribution of halos identified from one realization of mock light-cone data generated stacking TV-body outputs. The large 
filled circles denote halos with M halo [h- 1 M Q ]> 3 X 10 14 , the small filled circles are for 3 X 10 14 > M ha i [ft- 1 Mg]> 8 X 10 13 and the dots are for 
8 X 10 13 > MhaicJ/i - 1 Mq]> 1 X 10 13 . Only the halos within redshift interval between 0.1 and 0.7 are displayed. Upper-right: The gray map shows the 
lensing convergence map on the 512 2 pixels with the pixel size of 0.24 arcmin obtained from ray-tracing experiment through the same mass distribution as 
used in the upper-left panel. The circles show the positions of halos with u n { w > 4, the center of which is set at the position of each halo and the radius of 
which shows the projected virial radius (see §. l3.4l for details). Bottom-left: Same as the top-left panel but the convergence map is smoothed by the Gaussian 
window function with 6q = 1 arcmin. The edge regions of the field with 5 arcmin width are masked because of incomplete smoothing due to the lack of k 
data on outer region. The gray scale of this plot is same as that of the bottom-right panel. Bottom-right: Same as the bottom-left panel but the Gaussian noise 
is added before smoothing is applied. From comparison with the other plots, it is clear that high peaks due to massive halos are still identifiable. 



which implies that an extended mass distribution is likely associ- 
ated with eccentric mass distributions such as halos having sub- 
structures or on-going merger. We note that the PDF of axial ratio 
peaks at a/c ~ 0.5 and thus most of halos significantly deviate 
from the spherical symmetry. This broad distribution leads to sig- 
nificant scatter in the relation between the expected and observed 
lensing signals as seen below. 

Figure fTol shows the PDF of the shape parameter R defined by 
eg. J23t for massive halos with M > 3 X 1O 13 /i _1 M . A broad dis- 
tribution is evidently seen and we checked that the distribution does 



not largely change over the redshift range we considered (z < 0.8). 
It is important to note that the median of the distribution is less than 
unity (~ 0.92), thus more than half of the massive halos (specif- 
ically 55% for halos with M > 3 x 10 13 h~ 1 M Q irrespective of 
the redshift) have R < 1. This is a natural consequence of the ran- 
dom orientation of the halos' major axis; among the spatial three 
directions, it is more likely that the major axis of a halo is elon- 
gated along either of two directions perpendicular to the line-of- 
sight rather than along it. The skewed distribution leads to system- 
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atic bias in the relation between the expected and observed lensing 
signals (see fl5.2l and Appendix Iclfor details). 



3.5 Visual impressions 

Figure \H\ represents the lensing convergence maps and the celes- 
tial distribution of halos, which are taken from one of the 50 re- 
alizations. In the top-right panel, the gray scale shows « value on 
the 512 2 pixels with the pixel size of 0.24 arcmin. The bottom-left 
panel shows the same convergence map as in the top-right panel 
but after smoothed by the Gaussian filter function of 9g = 1 ar- 
cmin. Adding the Gaussian noise yields the noisy map shown in the 
bottom-right panel. While the noise indeed alters the appearance of 
the convergence map, the very high peaks remain identifiable. 

The top-left panel shows the distribution of all the halos iden- 
tified by the FOF scheme in the original iV-body simulations from 
which the convergence maps in the other plots are generated. Note 
that the size of point scales with halo masses. A good correspon- 
dence in positions between massive halos and high lensing signals 
on the k map is seen. This correspondence is more clearly seen 
if halos are selected by means of the expected J/ n f w value: in the 
convergence maps (the top-right and two bottom panels), the posi- 
tions of halos with ^ n f w > 4 are marked by the circles with their 
centers being placed at the halo center-of-mass position, and their 
radii show the angular virial radius (see i|3.4t . The improvement in 
the correspondence reflects the fact that the lensing signal is deter- 
mined by combined contributions of the halo mass, the redshift and 
the halo concentration, rather by the mass alone. 



4 RESULTS ON PEAK COUNTS 

In this section, we study the number counts of high peaks in the 
lensing convergence map. To do this, we use three kinds of k maps; 
(i) the noise-free k maps, (ii) the noisy k maps where the Gaus- 
sian noise was added in the manner described in £13.31 and (iii) the 
pure noise maps without the cosmological lensing signals. All these 
maps are smoothed with a Gaussian window function with 9a = 1 
arcmin (eq. 1111 ). We identify positive (negative) curvature peaks in 
the maps as pixels that have higher (lower) values of k than all the 
surrounding eight pixels (Jain & Van Waerbeke 2000, and also see 
Miyazaki et al 2002 for the first application to actual observational 
data). The peak height is defined as the ratio of the convergence 
value at the peak to the noise RMS, as given by eq. 1 141 . For all the 
= 0.02 (<7 e = 0.4 and n g — 30arcmin -2 ) is 



0.6 



three cases, a n 
adopted. 



4.1 The peak PDF 

The probability distribution function (PDF) of peaks in the con- 
vergence map is studied in detail by Jain & Van Waerbeke (2000). 
Here we briefly summarize the characteristics of the peak PDF that 
are relevant for the subsequent discussion. 

Figure fT2l shows the peak PDFs in the noise-free k (dashed 
histogram), the noisy k map (solid histogram) and the pure noise 
map (dotted histogram). The upper panel shows the sum of the pos- 
itive and negative curvature peaks, while the lower panel shows 
the contribution from the positive peaks only. In the absence of 
the noise, the peak PDF reflects the characteristics of the nonlin- 
ear gravitational clustering. Most of negative peaks are due to the 
negative curvature peaks (troughs). There is a cutoff in the negative 
peak height because the mass density contrast cannot be lower than 




Figure 12. The simulation results for the PDFs of the peaks in the noise- 
free k (dashed histogram), in the noisy re map (solid histogram) and in the 
pure noise map (dotted histogram). Shown is the cumulative results from 
the 50 realizations. The upper panel shows the sum of the positive and neg- 
ative curvature peaks, while the lower panel shows the contribution from 
the positive peaks only. In the total 50 X 2.048 2 square degrees' area, there 
are 59,908, 107,165 and 114,250 peaks for the noise-free re, the noisy re 
map and the pure noise map, respectively. The corresponding numbers for 
positive curvature peaks are 31,043, 53,123, and 57,141. 



— 1. One can also find a tail toward the high positive peak height 
which arises from the lensing signals due to massive halos (Kruse 
& Schneider 2000). Another noticeable feature seen from the lower 
panel is that some of the positive curvature peaks have a negative 
peak height, which occurs when the peak due to a halo is super- 
posed on a void region in the same line-of-sight. If the noise field is 
Gaussian, the PDF from the pure-noise map is derivable from the 
the Gaussian random field theory (Bardeen et al. 1986; Bond & Efs- 
tathiou 1987; Van Waerbeke 2000). The symmetric bimodal feature 
of the peak PDF is due to the superposition of two contributions 
from peaks with positive and negative curvatures. While adding the 
noise to the cosmological ft map alters the peak PDF drastically, 
there remain some characteristic features that contain useful cos- 
mological information. First, the asymmetric double peaks are due 
to the skewed feature of the cosmological re PDF as stated above. 
Second, there are excesses over the noise PDF on both high pos- 
itive (v > 3) and negative peaks (v < —3). These features were 
indeed detected in the recent weak lensing survey by Miyazaki et 
al. (2002). 

It is worth pointing out that the pure noise PDF can be directly 
estimated from an actual observed data. First, we can generate a 
pure noise map by applying mass-reconstruction to a randomized 
catalog where the orientation of source galaxy images is random- 
ized (to remove coherent image distortions from lensing signals). 
Then we can measure the peak PDF from the pure noise map, as 
was performed in Miyazaki et al. (2002). This allows one to esti- 
mate the contamination rate (the false peak rate) directly from the 
data. 
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Figure 13. Counts of peaks per 1 square degrees. The histograms are from 
the numerical experiments; the noisy k maps (solid), the noise-free k maps 
(dashed) and the pure noise realizations (dotted). The long-dashed line 
shows the model prediction, eq. 1161 with eq. 1151 . and the dot-dashed line 
shows the prediction by Gaussian random field theory (e.g., Bond & Efs- 
tathiou 1987; Van Waerbeke 2000). 



4.2 High peak counts 

In the following, we restrict our discussion to the number counts of 
high peaks that are mainly due to massive halos. Figure fl3l shows 
the number density of (positive curvature) peaks per 1 square de- 
grees. The histograms show the average values computed from the 
50 realizations of the numerical experiments; the noisy « maps 
(solid), the noise-free k maps (dashed) and the pure noise maps 
(dotted). The long-dashed line shows the model prediction from 
eq. 1 1 61 . while the dot-dashed line shows the prediction for the pure 
noise counts (Van Waerbeke 2000). From comparison of the results 
for the lensing maps with and without the noise, it is clear that the 
noise significantly increases the peak counts over a range of peak 
heights we have considered. Even high peaks with v > 5 are simi- 
larly boosted, while there is no such event in the pure noise counts, 
as shown by the dotted curve. There are two main reasons for the 
enhancement (see £|5,3l and Appendix [B] for further details). First, 
addition of the noise leads to modifications of the peak height as 
well as of the peak position. Second, the noise generates spurious 
peaks in the noisy maps. In particular, the latter effect is likely to 
show up in a biased way depending on the peak height. Since the 
convergence field around a given high peak has greater amplitude 
than the average field, adding the noise onto such a region is more 
likely to induce spurious high peaks around the real peak. We de- 
velop a correction scheme to the model peak counts in Appendix 
IrTl The corrected counts are in reasonable agreement with the sim- 
ulation result for high peaks v > 5. 



5 RESULTS ON THE CORRESPONDENCE BETWEEN 
HALOS AND PEAKS 

In this section, we examine correspondences between massive ha- 
los and peaks in the weak lensing convergence maps, using the re- 



Table 3. Classification of matching relations between halos and peaks. The 
double-lined arrows denote the primary match, while the single-lined ar- 
rows denote the secondary match. The last two columns summarize the 
number of halos/peaks in each class, the third (fourth) column is for match- 
ing between the halo catalog and the peak list from the noise-free (noisy) 
n map. Note that halos with i^ n f w > 3 (856 halos in total) and peaks with 
^lcns > 3 (669 peaks) and f no i S y > 3 (2417 peaks) are considered. 



class 


matching relation 


number of objects 






^nfw'^lcns 


^nfw'^noisy 


i 


halo^S peak 


475 


624 


ii 


halo with no paired peak 


367 


222 


iii 


peak with no paired halo 


177 


1637 


iv 


halojzf; peak 


14 


10 


V 


halo^Z! peak 


17 


156 



suits of our ray-tracing experiments and the halo catalogs created 
from the A r -body simulation outputs. In doing so, we follow a re- 
alistic procedure which is to be used in weak lensing surveys: we 
focus on high peaks above a certain threshold in the lensing map as 
strong candidates for halos. The usefulness of the weak lensing halo 
search can be quantified in terms of completeness and efficiency. 
As will be shown later in detail, choosing a threshold v t h — 4 — 5 
gives an optimal balance between the efficiency and completeness. 
We will also carefully examine if there is a correlation between the 
detectability of halos and their shape, because it could produce a 
bias in halo catalogs constructed from the weak lensing method. 

In order to make realistic weak lensing survey catalogs, we 
impose a criterion on the field selection of simulations so that no 
nearby massive halo is included in the target fields. Specifically, 
we discard simulated fields if there is a halo (or possibly more) 
with A/haio > 1O 14 /i -1 M0 at z < 0.1. Similar criterion, such as 
no known nearby X-ray luminous cluster of galaxies or no nearby 
optically identified massive cluster, is usually imposed on the field 
selection of an actual observation, because such a nearby massive 
cluster occupies a large portion of a survey field and makes a distant 
halo search difficult and inefficient. 3 1 out of the 50 simulated fields 
pass this criterion. In addition, we discard one field in which there 
is a massive halo with v n tw > 5 having flag = 2 (see £13.41 for 
a description of flag) 3 . We analyze the remaining 30 fields. We 
note that we do not use data within 5 arcmin of the edge of each 
simulated field because of the incomplete smoothing due to the lack 
of data at the outer region. The total field size is then 106 square 
degrees. 



5.1 Matching halos with peaks in lensing mass maps 

We shall now examine relations between the peak height in weak 
lensing maps and the properties of halos through the expected peak 
height for the NFW halos (^ n f w , see H2.lt . We restrict our discus- 
sion to halos with i^ n f w > 3 only. The halo catalogs contain 856 
such halos in total. On the other hand, the ray-tracing simulations 
directly enable us to find peaks in the n maps. We consider two 
kinds of the k maps; those with and without the noise due to the 
intrinsic ellipticities. The peaks in the two maps are identified in 
the manner described in S|4| The noise-free maps contain 669 peaks 
with u\ cnB > 3, while the noisy maps contain 2417 peaks with 



3 Including this field in our analysis below does not change the conclu- 
sions, but we decided to discard it for having conclusive results. 



© 2003 RAS, MNRAS 000.[T1|22| 



12 T. Hamana, M. Takada & N. Yoshida 



^noisy > 3. Thus, the noise significantly enhances the peak counts. 
It is also worth noting that the 669 peaks in the noise-free maps is 
significantly smaller than the expected number of peaks from the 
halo catalogs. We will carefully investigate this discrepancy in the 
subsequent sections. 

We consider halos and peaks with v > 3 only in order to 
avoid mismatches between peaks and physically unrelated halos. 
Because of large number density of low mass halos, including such 
halos rapidly increases the chance of mis-matching. This is espe- 
cially the case when one attempts to match peaks in the noisy k 
maps with halos because false peaks due to the noise can be ac- 
cidentally matched with physically unrelated halos. The criterion 
v > 3 is chosen based on the fact that most of such high peaks are 
expected to be caused by real halos, even when the lensing map is 
contaminated by the noise (see Figure [T3i . Also the criterion value 
is reasonably lower than the observational threshold value which is 
generally set by v > 4. This guarantees most halos of interest to us 
are included in the final halo catalog. 

We carry out the peak-halo matching in two directions: 
whether a peak has a corresponding halo in the halo catalog, and 
whether a halo has a counter-peak in the simulated « map. Note 
that the peak position does not necessarily agree with the halo's 
center-of-mass position exactly and this is more likely the case for 
halos that have substructures or are in the merger process. For this 
reason, we search for a matched pair candidate within a radius of 
12 pixels (2.88 arcmin) from the peak position or from the halo 
center. This maximum angular separation 2.88 arcmin is chosen so 
that it is larger than the smoothing radius of 9q = 1 arcmin, while 
it is still smaller than the angular virial radius of a massive halo at 
0-1 z ^ 0.7. If there are more than one pair candidate within 
the radius, we take the closest one as the primary candidate and 
consider the others as the secondary candidate. 

We classify the peak-halo matching into the following five 
classes: Class-(i) is that the both directions agree; both a halo and 
a peak are their primary pair candidate. Class-(ii) is defined as ha- 
los that have no paired peak within a radius of 12 pixels from the 
halo center. We will call such halos as "missing halos". Class-(iii) 
is defined as peaks with no paired halo within a radius of 12 pixels 
from the peak position. We will often refer to such peaks as "false 
peaks". Class-(iv) denotes a halo that has a matched peak candi- 
date, but is not the primary matched halo of that peak. Class-(v), 
which is the opposite to the class-(iv), is that a peak has a matched 
halo candidate but is not the primary matched peak of that halo. 
The resulting numbers for the noise-free and noisy maps are sum- 
marized in Table|3| 

Typical properties of halos or peaks in the class-(ii-v) are ex- 
plained as follows. The class-(ii) halos are those which have a peak 
height lower than the criterion value (v — 3) though the expected 
peak height, f n f w , exceeds it. This is mainly due to departures from 
the NFW profile or/and the projection effect. In the class-(iii) peaks 
there are two kinds: The first is the peaks which have corresponding 
halos but have the expected peak heights, v^s™, lower than the cri- 
terion value, though the measured peak heights in a k map exceed 
it. This is again due to departures from the NFW profile or/and the 
projection effect. The other is the peaks which do not have corre- 
sponding halos but are generated by the noise. The former accounts 
for most of the class-(iii) peaks in noise-free k maps, while the lat- 
ter accounts for a large fraction of those in noisy k maps. The class- 
(iv) halos are due to the projection effect; if two or more halos are 
located at different redshifts along almost the same line-of-sight, 
after smoothing of the n map only one peak appears, because the 
smoothing merges those into one structure. Finally, in the class-(v) 
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Figure 14. The correspondences between halos and peaks identified in the 
noise-free re maps in I'nfw-^lens plane. The class-(i-iii) objects are shown. 
The class-(ii) objects are placed within a narrow track around v\ BUS ~ 2, 
regardless of their true value, and the class-(iii) objects are placed at u n { w ~ 
2 as well. 

peaks there are two kinds: The first is the peaks due to substruc- 
tures within a halo. This accounts for most of the class-(v) peaks 
in noise-free k maps. The other is false peaks due to the noise gen- 
erated around the primary peak due the massive halo. Because the 
convergence field around the high peak is biased and thus larger 
than the average field, adding the noise into such region is more 
likely to generate false peaks with a relatively high peak height. 
This accounts for most of the class-(v) peaks in noisy k maps. 

In the rest of this section, we consider the class-(i-iii) objects 
only, and do not deal with the class-(iv) and (v) objects, because 
the class-(iv) and (v) objects are the secondary matched pairs asso- 
ciated with the primary matched pair in the vicinity. This roughly 
mimics the actual observation procedure in which the secondary 
objects are investigated in detail in course of a follow-up observa- 
tion of the primary matched object. 

5.2 Halos versus peaks in noise-free k maps 

Let us first look into the correspondence between high peaks identi- 
fied in the noise-free k maps (lAens) an d properties of halos through 
the expected peak height (^ n f w ). The main purpose of this subsec- 
tion is to examine the amplitude of the irreducible scatter in the 
relation induced by both the individuality of halo's mass distribu- 
tion and the projection effect. 

Figure 1141 shows the distribution of the class-(i-iii) objects 
in the f n fw-ficns plane. Since the class-(ii) halos do not have a 
matched peak with height Ui ens > 3 in the k map (although the 
halos have a high expected peak height of i/ n fw > 3), such halos 
are placed in a narrow locus around v\ ens = 2 for illustrative pur- 
pose. For the similar reason, the class-(iii) peaks are placed around 
f n fw = 2. The distribution of the class-(i) halos shows a large scat- 
ter around the ideal relation ^i ons = fnfw (shown by the solid line 
in the upper panel) over the range of v we consider. It is worth 
noting that the distribution for a fixed i/ n f w value is slightly biased 
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Figure 15. The scatter plot for the distribution of the halos in z-Q / Qnpw 
plane (lower panel) and in z—R plane (upper panel). The filled circle and 
triangle symbols are for the class-(ii) halos with ^ n f w > 5. For comparison, 
the class-(i) halos with 5 < !^ n fw < 7.5 are also plotted by the cross 
symbols. It is apparent from the lower panel that the class-(ii) halos can be 
separated into two subclasses: those with Q/Qnpw > 2 are plotted by the 
filled triangles, while the rests are plotted by the filled circles. 



Figure 16. Upper Panels: The completeness f c (y > fthl^nfw > v min) 
defined by eq. 1241 . as a function of the threshold peak height u t h- Three 
cases for the minimum expected peak value, u m i D = 3, 4 and 5 are shown 
by the filled squares, triangles and circles, respectively. Lower Panels: The 
efficiency /effnfw > u min\ u > u th) defined by eq. 1251 . as a function 
of the threshold peak height u t ^. The left panels are for the noise-free case 



(u = i^i cns ), while the right panels are for the noisy case (u - 



toward lower v\ ens values. The mean and RMS of the differences, 
Vlens — iAifw, are —0.24 and 1.2, respectively. Since no noise is 
added to the k map in this case, the scatter and bias must be as- 
cribed to the two effects; (1) the mass distribution within a halo is 
not spherical nor completely universal, causing deviations of the 
peak heights from the expected values computed from our simple 
model, and (2) a chance projection of structures along the same 
line-of-sight alters the lensing signal due to the halo. For the for- 
mer, although we have so far used the average relation for the halo 
concentration, given by eq. J2j, it was shown in Jing (2000) that 
the halo concentrations among halos with a given mass obey the 
log-normal distribution with the variance cr(lnc n f w ) m 0.2. This 
roughly corresponds to a 20% RMS scatter in the c n f w relation. 
Given the approximate relation k oc c n ~[ w ' for the peak height due 
to a massive halo (see J2j, the RMS scatter in the c n fw distribution 
translates into 20 — 30% RMS scatters in the lensing signal i/ n f w 
and therefore equivalently in the ^i cns -^nfw relation. For example, 
for halos with ^ n f w = 5, the RMS of <j„ ~ 1 is caused. Sec- 
ond, the lensing projection due to physically unrelated structures 
along the same line-of-sight sight of a halo also contributes to the 
scatter in ^i ens - Following Hoekstra (2001), this scatter can be esti- 
mated from the cosmic weak lensing variance under the assumption 
that the halo and the large-scale structures are uncorrelated. For the 
Gaussian smoothing scale 9q = l', the amplitude of the (linear) 
convergence variance is <r K (l') = 0.01 for the ACDM model, cor- 
responding to the RMS deviation of a v =0.5 because of the noise 
0"noise(l') = 0.02 (see eq. 1121 ). Hence, these two effects can basi- 
cally explain the amplitude of the scatter in the Hens-t'nfw relation 
in Figure [141 

We also argue that the systematic bias in the v\ eiis — i^ n f w rela- 
tion can be explained by the following two effects. The first is the 
fact that more than half of the halos have R < 1 as discussed in 



i|3.4l and Figure [Tol Since R positively correlates with i^i ene /vnfw 
(see Figure IC2l in Appendix Icl. it causes the bias. The second is 
the projection with a void region in the same line-of-sight which re- 
duces the lensing signal (see Metzler et al. 2001; White et al. 2002; 
Padmanabhan et al 2003 for detailed study on the projection effect). 

Let us turn to the class-(ii) halos in Figure fl4l Among the 
367 class-(ii) halos, there are 17 halos that are supposed to produce 
peaks with v > 5 in the k maps if the halos have the NFW pro- 
file (i.e., !/ n f w > 5 in our notation). Why are these halos missing? 
Figure \^\ explains this in terms of deviation of the mass distribu- 
tion within a halo from the spherically symmetric NFW profile. 
We again use the parameters Q and R to quantify the mass dis- 
tribution within a halo (see equations 1201 and 1231 . respectively): 
Q/Qnfw < 1 (> 1) characterizes a halo that has centrally more 
(less) concentrated mass distribution than expected from the NFW 
profile, while R > 1 (< 1) describes that the mass distribution is 
elongated along (perpendicular) to the line-of-sight (R — 1 for the 
spherical symmetric halo). Therefore, a halo with Q/Qnfw > 1 
and/or R < 1 likely leads to a smaller lensing signal than expected 
from the NFW profile (see Appendix |C| for more detailed analy- 
sis of the relations between the halo shape and the amplitude of 
the lensing signal). The upper panel of Figure fT5l shows the dis- 
tribution of the class-(ii) halos in z-Q/Qnfw plane, while the 
lower panel shows the distribution in z-R plane. For comparison, 
the figure also shows the distributions for the class-(i) halos with 
5 < f n fw < 7.5 by the cross symbols. From the lower panel, it is 
clear that the class-(ii) halos can be divided into two distinct sub- 
groups; those which have Q/Qnpw > 2 (plotted by the filled tri- 
angles) or Q/Qnfw < 2 (plotted by the filled circles). We indeed 
find that halos with Q/Qnfw > 2 show clear signatures of sub- 
structures or on-going merging, resulting in the extended mass dis- 
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Figure 17. The scatter plot for the halos in the redshift-halo mass plane. The contours show the expected peak height ^ n f w from 3 to 9 at an interval of 
At/ n f w = 1. The class-(i) halos matched with peaks identified in noise-free k maps are plotted in the upper-left (v\ cns > 5), upper-right (5 > v\ cus > 4) and 
lower-left panels (4 > !^ cns > 3). The lower-right panel shows the scatter plot for the class-(ii) "missing halos". 



tribution (the k map of such halos is shown in Figure IaTI for exam- 
ple). The remaining objects with Q/Qnfw < 2 tend to be located 
at a low redshift (z < 0.25), where the lensing efficiency is rela- 
tively small for the source redshift z s — 1 we consider, and have 
a small R, mostly R < 1. The examples are shown in Figure s lA2l 
and lA3l These halos typically have prominent substructures and the 
whole mass distribution within the halo is indeed elongated along 
the perpendicular direction to the line-of-sight. Except for such ir- 
regular systems, almost all massive halos with f n fw > 5 generate a 
high peak with z^i cns > 3. It is important to note that the distribu- 
tion of the class-(ii) halos in Figure fT5l is not clearly distinguished 
from that of the class-(i) halos, reflecting the fact that the lensing 
signal is also affected by the projection effect. 

Turn next to the class-(iii) objects in Figure 1141 which are 
the peaks with ^i cns ^ 3 in the lensing map, but do not have 
the corresponding halos in the halo catalog (contains only halos 
with z/ n f w > 3). Among the class-(iii) objects there are 8 events 
with fi cns > 5 which we call "false peaks". We carefully looked 
into each case in the k map and searched for a possible related 
halo in the halo catalog. We found that the 5 events are associated 
with substructure of massive halos or the secondary local peaks 
of merger systems (note that the peaks are separated from the pri- 
mary peak of the halo center by more than our searching distance, 
12pixels = 2.88 arcmin). Figure |A4| shows the k map of one of 
such cases. The remaining 3 peaks are due to a chance projection 
of two moderately massive halos at different redshifts (Figure lA5l 
shows the one example). Therefore, except for the three chance pro- 



jection systems, the false peaks with i^ ens > 5 are all associated 
with a real massive halo as the secondary peak. 

In order to quantify the ability of weak lensing survey to locate 
massive halos, we consider the completeness and efficiency. We de- 
fine the completeness by the fraction of halos, which are detected 
as high peaks above a given threshold vth in the k map, relative to 
all the halos with the expected lensing signal, t/ n f w , greater than a 
given minimum value f m i n : 

f , | . s Ni{V\ ens > t-thKfw > JAnin) 

Jc{lAens > fth|fnfw > Vjoia) = — r , (24) 

JVi+ii(^ n f w > V m \ 

where JVi+ii(i/ n f w > i/ m i n ) denotes the number of halos with 
^nfw > fmtn which, in our classification, corresponds to the num- 
ber of the class-(i) and (ii) such halos, and Ni(y\ CTts > fthl^afw > 
;y m i n ) denotes the halos having ^ n f w > v mlrL matched with the 
peaks with i4 e ns > v th (thus these are the class-(i) objects). Sim- 
ilarly, we define the efficiency by the fraction of peaks, which are 
associated with halos with the expected peak height exceeding a 
given minimum value ^ m i n , relative to all the peaks with height 
greater than a given threshold f t h- 

, , \ Ni(u n f w > fminl^lcns > fth) 

M^nfw > fminl^lens > ^th) = T7 } ; , (25) 

JVi + iii(fl ons > U th ) 

where A r i + iii(i/i ens > f t h) denotes the number of peaks with 
!4ens > Vtb which corresponds to the number of the class-(i) and 
(iii) such peaks. 

The results are plotted in the left panels of Figure ll6l as a func- 
tion of f t h- The upper panels show the completeness (/ c ) and the 
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Figure 18. The correspondences between the peak heights identified in the 
noise-free and noisy k maps. The peaks, which have v\ CYls > 3 but are 
missed in the noisy n map or have v nolsy < 3, are placed within a narrow 
locus at Vnoisy ~ 2, regardless of their true u lYolay value. Similarly, the 
peaks, which have v ao \ sy > 3 but no paired peak with i^i cns > 3, are 
placed at ui cns ~ 2, 

lower panels show the efficiency (/ e ). As expected, these two es- 
timates behave inversely, that is, when f E goes up, f c goes down. 
Therefore the optimal choice of the i/ t h value should be determined 
by balance between the two estimates, which we discuss in the lat- 
ter section with taking into account effects of the noise. 

The plot for the completeness shows that only about half of 
halos with v n f w > 5 produce peaks with !/i cns > 5 in the noise- 
free k maps, and that 77 percent of the halos have viens > 4. The 
completeness becomes worse if Vth. > ^min, as expected. It is also 
shown that a good completeness (say f c > 0.8) is attainable if one 
sets f t h ;$ ^min — 1. This can be understood from the fact that the 
RMS scatter caused by the individuality of halos and the projec- 
tion effect is a v ~ 1.2. Note that some fraction of the halos with 
i/iens much smaller than the expected peak heights f n fw are irregu- 
lar systems or un-relaxed systems, and therefore the completeness 
is likely to be improved for samples of regular, well relaxed halos. 
It should be also noted that the noise increases the completeness 
(see il5.4l and the upper right panel). 

Let us turn to the efficiency. Since no noise is added to the k 
maps, the efficiency given here should be regarded as an estimate 
for an ideal lensing survey. Clearly, the efficiency is very high. If 
one sets the detection threshold to be — 5, 73 percent of all 
the peaks correspond to real halos with «/ n f w > 5, and 95 percent 
for i/ n f w > 3. Even for v t h = 4 the efficiency is still high, 90% 
for i/ n f w > 3. Therefore, in the absence of the noise, the contami- 
nation rate (the fraction of false peaks) is very small for detection 
threshold i/ t h > 4. 

Figure fTTl shows the scatter plot of the class-(i) and -(ii) halos 
in the redshift-halo mass plane. The contours show tracks of the 
expected peak height ^ n f w from 3 to 9 at an interval of Af n f w = 1. 
The class-(i) halos matched with peaks in noise-free k maps are 
plotted in the upper-left (^i cns > 5), upper-right (5 > V\ e ns > 4) 
and lower-left panels (4 > Viens > 3), respectively. The lower-right 
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Figure 19. Same as Figure fHl but for the correspondence between halos 
and peaks identified in the noisy k maps. 

panel shows the scatter plot for the class-(ii) halos, namely "missing 
halos". Importantly, halos within a small range of the peak height 
are spread over a wide range of masses due to the projection effect 
and the individuality of the halo mass distribution. 

5.3 Effects of the noise on the lensing peaks 

Before examining correspondences between peaks identified in the 
noisy k maps and halos in the halo catalog, we quantitatively exam- 
ine how the noise alters the peaks due to halos and generates false 
peaks in the lensing map. 

Matching between the noise-free peaks and the noisy peaks 
was done in the same manner described in 85.11 In what follows, 
the peaks classified in the class-(i-iii) (its analogy to ficns-^noisy 
correspondence) are considered. Again, we only deal with peaks 
with v > 3. 

Figure [T8l shows the relations between i^ ens and fnoisy plotted 
in the same manner as in Figure^] We emphasize that the statisti- 
cal properties of the fi cn s-fnoisy relation do not follow the Gaussian 
distribution, even though the Gaussian noise is added. The reason is 
as follows. The noise not only alters the peak height from massive 
halos but also changes the peak position. For example, if a large 
positive noise is accidentally added on a neighboring pixel around 
the original peak position, the pixel could be identified as a peak 
that is even higher than the original peak. Therefore, the distribu- 
tion of i^noisy's is not simply expressed as the convolution of that 
of f lens's with the Gaussian noise, but the ficns-fnoisy relation be- 
comes asymmetric and is biased toward large i/ n oisy values. This 
explains the asymmetric lAens-^noisy relation clearly shown in Fig- 
ure^] The mean and RMS of the differences, tWsy— lAens, among 
the class-(i) objects with fi ens > 4 (this restriction is imposed in 
order to avoid an incomplete sampling of low ^noisy pairs) are 0.31 
and 0.90, respectively. Therefore the peak heights of massive halos 
identified in a noisy k map are, on average, boosted by the noise. 
This bias partly accounts for the excess in the peak number counts 
over the simple theoretical prediction, eq. 1161 . as discussed in 84.21 



16 T. Hamana, M. Takada & N. Yoshida 




The class-(ii) objects here correspond to the missing peaks, the 
peaks that have viens > 3 but are absent in the corresponding noisy 
k map or have v uo i sy < 3. Except for one peak with ^ ons = 5.1, 
the missing peaks sit in the range of vu>na ^ 4. These peaks are 
likely to be erased or their heights are decreased by the addition 
of the negative noise. The one exception (the peak with ^i cns = 
5.1 at the missing peak locus) is a local peak associated with a 
substructure of a massive halo which disappears in the noisy k map. 
We may, therefore, safely conclude that more than 80% of the peaks 
with view ^ 3 in the noise-free map remain in the noisy k map. 

Next, we consider the false peaks, the peaks that are identi- 
fied in the noisy « map with ^noisy > 3 but does not correspond 
to peaks in the noise-free « map or have ui enB < 3, which are 
placed around vi cns ~ 2 in Figure [1*81 A large fraction of the false 
peaks, especially for ^noisy < 4, are the false peaks generated by 
the noise (see Sj4]for details), and the small fraction (especially 
those with relatively high ^ no isy) is associated with real peaks on 
which a large noise is accidentally added. Actually, there are 14 ob- 
jects with t'noisy > 5 in the/a/.se peaks locus. All of those have a 
corresponding peak in the noise-free n map with low peak height 
1.5 < ficns < 2.9 on which a relatively large noise of 2 ~ 3.5<r is 
accidentally added. Therefore, all peaks with ^ no isy > 5 identified 
in the noisy k map are associated with real peaks in the noise-free 
k map. On the other hand, most of the class-(iii) peaks with lower 
heights are false peaks generated by the noise. In fact, one can find 
the rapid increasing number of the false peaks toward low peak 
heights, especially for ^ no isy < 4, as shown in 21 It immediately 
follows from this that in order to reduce the contamination rate, a 
high threshold value, u t h ^ 4 at least, is required. 



5.4 Halos versus peaks in noisy k maps 

In this subsection, we examine correspondences between halos in 
the halo catalog and peaks in the k map including the noise due 
to the intrinsic ellipticities. Figure [T9l shows the distribution of the 
class-(i-iii) objects in v n fw-V noisy plane plotted in the same manner 
as in Figure 1141 The comparison with Figure [2| reveals that the 
noise not only increases the scatter in the fnoisy-i'nfw relation but 
also causes the systematic bias in that the peak heights in the noisy 
map are statistically boosted, as seen in the previous subsection. 
Interestingly, this bias leads the completeness for the weak lensing 
halo search to be improved. The upper-right and lower-right panels 
in Figure [T6l show the completeness and efficiency, defined in the 
same manner as eqs. I24i and 1251 . for the noisy maps. About 81 
(63) percent of massive halos with u n {„ > 5 (> 4) are identified as 
peaks with f no isy ^ 4 in the noisy maps. This is one of the most 
important results of this paper, indicating that a high completeness 
of the weak lensing halo search could be attained in a realistic data. 

There are 13 missing halos that have ^ n f w > 5 but v no isy < 3 
or have no corresponding peak in the noisy k map. We look into 
each case in the k map and find that 1 1 halos out of 13 are irreg- 
ular systems such as un-relaxed systems having substructures or 
on-going merger. Therefore, the completeness of the weak lensing 
halo search for regular, well relaxed systems would be higher than 
the estimate plotted in Figure [T6l The remaining 2 missing halos, 
which have v n fw = 5.3 and 5.4, are halos which accidentally meet 
a large negative noise so that their k values are below ^ no isy = 3. 
The fraction of such halos over all halos with ^ n f w > 5 is 0.8% (2 
out of 250 halos), which is indeed very small. Hence, we can safely 
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Figure 21. Same as Figure [l9l but the halo mass is used instead of the ex- 
pected peak height of v n f w . Note that only halos with u n f w > 3 are con- 
sidered. 



conclude that almost all halos with ^ n f w > 5 generate peaks with 
fnoisy > 3 even in the presence of the noise. 

Let us turn to the class-(iii) objects. Since the noise generates 
false peaks in the «: map, the efficiency becomes worse than the 
noise-free case (see Figure fT6l . The fraction of contamination by 
the false peaks is larger for a lower peak height because the noise 
generates relatively lower peaks, mostly u no i sy < 4. Therefore, if 
vth = 3 is employed as the detection threshold, only 28% of all 
peaks are associated with halos with j/ n f w > 3. Note that since we 
only consider halos with i/ n f w > 3, some fraction of missing peaks 
could be associated with halos with u n f w < 3. On the other hand, 
the efficiencies for higher detection thresholds are reasonably high: 
for fth = 5 (v t h = 4), 85% (56%) of all peaks with f no i sy > vth 
are associated with halos with j/ n f w > 3. 

From Figure [T3I it is clear that a better efficiency achieved 
by setting a higher detection threshold is a "trade-off" for a worse 
completeness. We also note that the peak counts varies sensitively 
with the threshold. For example, the counts N{y > 4) is about 
twice as large as N(i> > 5) (see §4). Hence, the detection thresh- 
old should be determined such that the efficiency and completeness 
are optimized for specific survey strategies. The results we have so 



far shown indicate that v no 



4 — 5 is reasonably optimal for 



currently typical weak-lensing surveys. To be more specific, the re- 
sults from the simulated k maps in the ACDM model indicate that, 
the mean number of halos with ^ n f w > 4 is A'haio(i / nfw > 4) = 37 
per 10 deg 2 , and on average 63% of such halos (23 halos out of 
37) are detected with ^noisy > 4. The completeness is better for 
higher i^ n f w halos: for halos with ^ n f w > 5, the mean number is 
iVhaioOnfw > 5) = 19 per 10 deg 2 , where 81% of the halos (15 
halos out of 19) are detected with /-'noisy > 4, and 12 halos have 
even high lensing signal of u no isy > 5. Concerning the efficiency, 
if one sets the detection threshold to be v t h = 4, the mean number 
of peaks is N vea k(v no \ ay > 4) = 62 per 10 deg 2 . Among the 62 
peaks, 23 are due to real halos with ^ n f w > 4, the 13 peaks are 
halos with 3 < u n { w < 4 and the remaining 26 peaks are likely to 
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Figure 22. The number fraction of halos identified by peaks with height 
higher than v tb relative to all the halos with mass greater than a given min- 
imum mass, M min , /m(tWsy > t'thl^halo > M min ) (see eq. |26|). 
Three cases of M min = 5 X 10 13 , 1 X 10 14 , and 2 X 10 14 /i~ 1 M Q are 
plotted. Note that only halos with u n f w > 3 are considered. 



be either halos with f n f w < 3 or false peaks due to the noise. On 
the other hand, if one adopts a higher detection threshold v t h = 5, 
the mean number of peaks is reduced to A r pea k(^ n oisy > 5) = 25 
per 10 deg 2 , however the completeness is improved significantly: 
in fact it is found that the 19 (3) peaks out of the 25 are signals from 
real halos with v n f v > 4 (3 < i^ n f w < 4). 

Figure l20l shows the scatter plot for the class-(i) and class- 
(ii) halos in the redshift-halo mass plane as in Figure fT71 From 
comparison with Figure fT71 one may find that more halos pass a 
certain detection threshold than in the noise-free map, because of 
the systematic bias caused by the noise. Another important point 
revealed in this figure is that the peak height in the noisy k map is a 
poor estimator of the halo mass. For example, masses of the halos 
detected with 4 < v no i S y < 5 span over an order of magnitude. 
Even in a narrow redshift range such as 0.3 < z < 0.4, there 
still exits a significant scatter in the halo mass. Note that we do not 
plot halos with v n f w < 3, which may cause a larger scatter toward 
small mass ranges. It should be also noted that in the distribution 
of the missing halos (lower-right panel) no strong tendency with 
redshift is seen, indicating that the completeness would not vary 
significantly with redshift. This is largely due to our definition of 
the completeness, that is, we define the completeness so that the 
effect of the selection function, which very strongly depends on the 
redshift, is decoupled. 

So far we have examined the relation between halos and peaks 
in a k map in terms of the peak height instead of the halo mass 
because the lensing observable is determined by combination of the 
halo mass, the redshift and the noise level, rather than by the halo 
mass alone. Nevertheless, it is interesting to see what mass range 
of halos produce the peaks in the noisy k maps, which is shown 
in Figure as in Figure 1191 Again only halos with is n f w > 3 
are considered. As can be seen, the peaks with a given height are 
caused by halos with a wide range of masses. We define the number 
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fraction on the analogy of the completeness defined by eq. 1241 but 
replacing the minimum peak height with the minimum halo mass: 

fm (^noisy > V t h | A/halo > Afmin) = 

iVi(^noi S y > ^thjA/halo > Mmjn) 
iVi+ii(M halo > M min ) 

where iVi+u(Mhaio > Af m i n ) is the number of (class-(i) and (ii)) 
halos with mass greater than M min , and Ni(v noiBy > f t h|Afhaio > 
Afmin) is the number of halos with mass greater than Afmin hav- 
ing peak height greater than v^. Figure l22l shows this estimate as 
a function of u^. It is apparent that f m declines toward higher 
detection threshold and for lower minimum mass. For detection 
threshold i/ t h = 4, more than 70 percent of halos with Afhaio > 
2 x 10 14 /i _1 Mq are identified, while about half of halos with 
Afhaio > 5 x W 13 h~ 1 M Q is identified, 



6 SUMMARY AND DISCUSSION 

We have investigated various aspects of the weak lensing cluster 
surveys employing both analytic descriptions of dark matter ha- 
los and the mock data of weak lensing surveys generated from 
numerical simulations. For the latter, we combined weak lensing 
ray-tracing through the mass distribution and the dark matter halo 
catalogs. Our major findings are summarized as follows. 

(1) In [J2| we examined the expected properties of weak lens- 
ing halos using the analytic descriptions of the dark matter halos, 
including the universal density profile (Navarro et al. 1996; 1997) 
and the Press-Schechter halo mass function (Press & Schechter 
1974: Sheth & Tormen 1999). We found that the Gaussian smooth- 
ing with 8q — 1 arcmin gives the largest expected weak lensing 
halo counts (which may however depend on the source redshift and 
noise properties). We computed the selection function of the weak 
lensing cluster survey and examined, in detail, the mass and redshift 
distribution of the weak lensing halos. It was shown that the de- 
tectability of halos depends not only on the mass but also strongly 
on the redshift. 

(2) In £13.41 we compared the model prediction of the weak 
lensing halo counts developed by Kruse & Schneider (2000) and 
Bartelmann et al. (2001) with the halo counts in the mock catalog, 
and found a good agreement. We also found a large scatter in the 
numbers of weak lensing halos within 4 square degree field, which 
is larger than the Poisson fluctuation. This can be explained by the 
strong clustering of massive halos. 

(3) In we tested the model prediction of the peak counts 
against the mock weak lensing data. It was found that a systematic 
bias is induced by the noise due to intrinsic galaxy ellipticities. This 
bias increases the peak counts, and hence the model prediction un- 
derestimates the counts. We developed a correction scheme in an 
empirical manner in Appendix [B] and showed that the improved 
model reproduces the peak counts reasonably well. 

(4) In EI5.2I using the mock weak lensing data combined with 
the halo catalog, we examined the matching between halos and 
peaks identified in the noise-free n map. This was done to clar- 
ify influence of the individuality of halos and the projection effect 
on a weak lensing halo search. We showed that these effects cause 
not only a large scatter in the ^i cn s-^nfw relation but also a sys- 
tematic bias. The mean and RMS of the differences v\ cnB — ^ n f w 
are —0.24 and 1.2, respectively. The level of scatter caused by the 
effects is comparable to that due to the noise. We argue that the 
negative mean value are due to a population of un-relaxed, less 
centrally concentrated halos and due to the fact that the more than 



half of the halos are elongated to the direction perpendicular to the 
line-of-sight. Also the projection of large under-dense regions may 
partially account for the scatter. It was also shown that the chance 
projection of massive halos in the same line-of-sight is very rare. 

(5) In £15.31 we used the noise-free k map and the noisy k, map 
to clarify how the peak distribution is affected by the noise and 
how false peaks are generated. We showed that the ^lcns-^noisy dis- 
tribution is biased toward large f no i S y values. The mean and RMS 
of the differences, ^ no isy — ^lons, among the matched peaks with 
lAenB > 4 and ^noisy > 4 are 0.31 and 0.90, respectively. Thus the 
noise not only generates the scatter but also systematically 'boosts' 
the peak heights. We found that almost all the peaks identified in 
the noise-free k map with Ui ens > 4 are identified in the noisy k 
map, thus the noise rarely erases the high peaks. All the peaks with 
^noisy > 5 identified in the noisy k map are associated with real 
peaks in the noise-free n map, whereas most of the false peaks (due 
to the noise) have a relatively lower peak height ^ I10 i S y < 4. 

(6) In £15.41 we examined the correspondence between the ha- 
los and peaks identified in the noisy k map. In particular, we studied 
the efficiency and completeness of the weak lensing halo survey. 
We found that the detection threshold of Vth — 4 — 5 gives an 
optimal balance between the efficiency and completeness. It was 
shown that about 81 (62) percent of massive halos with ^ n f w > 5 
(> 4) are identified as high peaks with vth = 4. This suggests that 
the completeness of the weak lensing cluster survey is reasonably 
high. Concerning the efficiency, we found that for i/th = 4 (5), 58% 
(82%) of all the peaks are real signals from halos with i/ n f w > 3, 
while 37% (63%) are signals from halos with higher expected peak 
height of i/ n f v > 4. Therefore, it is possible to attain relatively 
high efficiency and high completeness by selecting a moderately 
high detection threshold. 

(7) In i|5.4l we conclude that, for detection threshold v t h = 4, 
about a half of halos with Afhaio > 5 x 10 13 h -1 Mq and more 
than 70 percent of halos with Af na i > 2 x 10 14 /i~ 1 Af© indeed 
produce the lensing signals with ^noisy ^ 4. 

Weak lensing cluster search technique explored in this pa- 
per can be directly applied to data obtained by on-going/future 
wide field surveys. Extending the pilot 2.1 square degree sur- 
vey (Miyazaki et al. 2002), a wide field weak lensing survey 
(Suprime33, PI: S. Miyazaki) is being carried out with the wide 
field prime focus camera on Subaru telescope, Suprime-Cam 
(Miyazaki et al. in preparation). The Suprime33 will cover 33 
square degrees in total with the limiting magnitude R = 25.5 
(which provides n g > 30 and (z 3 ) ~ 1). For the survey exposure 
time, 1800 sec, and the field-of-view of the Suprime-Cam, 0.25 
square degrees, the effective survey cost is 2.5 hours per 1 square 
degrees including the overhead. The expected number of clusters to 
be located by the Suprime33 is 70 (120) for the threshold u& = 5 
(vth — 4). Therefore weak lensing surveys exploiting a wide field 
camera on a large telescope offer a reasonably efficient way to lo- 
cate massive clusters. The CFHT Legacy Survey 4 will significantly 
enlarge the survey area. Its "wide" survey will observe 170 square 
degrees in total with the limiting magnitude i' — 25.5 (which pro- 
vides n g ~ 30 and (z s ) ~ 1). The expected number of clusters is 
360 (600) for u t h = 5 (vth = 4), and thus it will provide invaluable 
data for studies on large-scale structure of the universe. Although 
we have primarily considered weak lensing surveys that are feasi- 
ble with current ground-based telescopes, future wide field surveys 
based on a space telescope will enable to accurately measure the 



4 http://www.cfht.hawaii.edu/Science/CFHLS/ 
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shape of distant very small galaxy images, which significantly im- 
proves the ability of weak lensing cluster surveys and allows to 
detect lower mass and/or higher redshift clusters. 

Finally, we note that it is in principle possible to enhance the 
S/N of the detection by choosing a suitable smoothing function 
(White et al. 2002; Padmanabhan et al. 2003). Padmanabhan et 
al. (2003) proposed the function of W(9) = (1 + 9/9 c )' 2 with 
9 C ~ 1 arcmin, which is motivated by the asymptotic behavior of 
the outer part of the projected NFW profile. Such optimized filters 
may improve the S/N if the matter distribution is indeed close to the 
NFW profile. On the other hand, its applicability to clusters whose 
mass distribution deviates from the NFW profile remains unclear. It 
needs to be explored how much improvement can be obtained both 
in efficiency and completeness. Optimizing the smoothing function 
clearly warrants further studies. 
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APPENDIX A: VISUAL APPEARANCES OF SYSTEMS 
WITH IRREGULAR SIGNALS 

Here we present some k maps, in which missing halos (the class- 
(ii)) or false peaks (the class-(iii)) with high signals {v > 5) exists, 
for examples of irregular systems. Fi gures |A 1 1 A5 1 show the (noise- 
free) convergence maps smoothed with Gaussian window function 
with 9a = 1 arcmin. Pluses denote the center-of-mass position of 
halos concerned and the circles show their virial radius. Crosses 
denote the position of peaks. 

Figures IA1IA3I show three examples of missing halos with 
the expected peak value ^ n f w > 5. Such high ^ n f w missing halos 
mostly have an un-relaxed appearance and are probably on-going 
merger systems. Figures lA4l and lA5l show two examples of false 
peaks with ^i cns > 5. Note that false peaks are rare in noise-free k 
maps. Most of false peaks with ^i cns > 5 are associated with sub- 
structures of massive halos. Occurrence of the chance projections 
is very rare, we found three cases with ^i cns > 5 in our 30 maps. 



APPENDIX B: SIMPLE CORRECTION SCHEME TO THE 
PEAK COUNTS 

As shown in fj4]the theoretical model of the peak counts given by 
eq. <16t underestimates the counts obtained from the realistic noisy 
k map by about 20 percent. We here discuss sources of this change, 
and develop the correction scheme to the peak counts in an empiri- 
cal manner. Here we focus on the high peaks of v > 5. 
The sources that modify the peak counts are: 

(a) The biased relation between i^nfw-^icns- In A5.2I it is found 
that the mean of the differences v n f w — v\ cns over class-(i) objects 
is —0.24. This means that the peak height computed assuming the 
NFW profile tends to overestimate the real peak height. The reason 
for this is discussed in A5.2I 



© 2003 RAS, MNRAS OOO.mSI 



20 



T. Hamana, M. Takada & N. Yoshida 



n 



*->? 

a ma 



i 



[arcmin] 



Figure Al. Missing halo, an example of a halo with extended matter dis- 
tribution (possibly in the process of merger). The plus denotes the center- 
of-mass position of the halo, and the circle shows the virial radius. The 
map apparently shows very eccentric density distribution. The halo has 
■Mhalo = 3.4 X 1O 14 /i _1 M at 2 = 0.37, and the shape parameters 
are Q/Qnfw = 2.1, R = 1.4, and c/a = 0.23. The expected peak value 
is ^ n f w = 7.0, but the v value measured in the lensing map is less than 3. 




Figure A2. Missing halo, an example of low- z un-relaxed system. The halo 
(Mhaio = 2.0 X 10 14 fo — 1 Mq, 2 = 0.13 and the shape parameters are 
Q/Qnfw = 1-2, R = 0.36, and c/a = 0.41) seems to be not yet reached 
to a relaxed stage, but probably is a on-going merger system. The expected 
peak value is u n { w = 6.6 but the actual v value is less than 3. 



(b) Spurious peaks generated by the noise. It is found in the fJ5| 
that spurious peaks yield about 10% fraction of the total number of 
high peaks (v > 5). 



(c) The scatter in the relation between v n 



y-V\ c 



due to the 



noise. As will be shown, this scatter increases the peak counts. 

(d) The biased relation between Vnoisy-iAew In 85.31 it is found 
that the noise produces not only the scatter but also the systematic 



bias, and that the mean of the differences u n 



: over class- 
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Figure A3. Missing halo, an example of an irregular system. The map ap- 
parently shows very eccentric density distribution. It is very likely that the 
halo (M halo = 1.9 X lO 14 /i _1 M , z = 0.20 and the shape parameters 
are Q/Qnpw = 1-3, R = 0.25, and c/a = 0.36.) is a on-going merger 
system. The expected peak value is u n { w = 6.1 but the actual v value is 
less than 3. 
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(i) objects is 0.3 1 . Therefore the peak heights in the noisy k map are 



Figure A4. An example of a local peak associated with the substructure of a 
massive halo. The cross denotes the position of the peak with fi ens = 5.3. 
This is the local peak associated with the substructure of the massive halo 
(A/ ha i = 3.5 x lO 14 h _1 M , z = 0.20, u nfw = 8.8 and the shape 
parameters are Q/Qnfw = 0.94, R = 0.61, and c/a = 0.66). The plus 
denotes the center-of-mass position of the halo, and the circle shows the 
virial radius. 



statistically boosted by the noise. The reason for this is discussed 
in g531 

First, we develop the correction scheme for the above point 
(c). To do this, we take a simplified way, employing the assumption 
that very high peaks are neither disappeared nor generated by the 
noise but their peak heights are altered by addition of the noise. In 
other words, we ignore, for the moment, spurious peaks due to the 
noise and the other sources. Under the above assumption, the peak 
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Figure A5. A very rare example of a chance projection of two moderately 
massive halos. The cross denotes the position of the peak with v\ cris = 5.1. 
This peak is not generated by a single halo but generated by the chance pro- 
jection of two moderately massive halos: (M halo = 7.5 X 10 13 ft _1 Af@, 
z = 0.35) and (M halo = 5.4 x 10 13 h- 1 M Q , z = 0.68). None of two 
halos alone can generate the lensing signal larger than v > 3. 
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Figure Bl. The number count of the lensing peaks per 1 squire degrees. 
The histograms are from the numerical experiments from the noisy k maps 
(solid), and the noise-free k maps (dashed). The long-dashed line shows the 
un-corrected model prediction computed from eq. 1161 . while the solid line 
shows the corrected model prediction (see text for detail). 



counts in the noisy k map is given by the convolution of the noise 
free PDF <16t with the Gaussian noise PDF: 



,{v)= I du' 



1 / v ri 

2 . exp r~ 



n pea k(^ - V). 



(Bl) 



We further assume that within a small range of v the peak counts 
in the original lensing map can be approximated by the exponential 



form np Ca k(V) = n„ exp(pz/), with a constant exponential index 
of p. As can be seen in Figure lBTl this would be a reasonable ap- 
proximation over a range of 3 < v < 8. This approximation allows 
us to compute the integration in eq. JB II analytically: 



dv' 



1 



2tt 



exp 



= exp (p 2 /2) n pcak (z. 



n* exp (p(v — u'y\ 
(B2) 



Thus the scatter due to the noise does increase the counts by a factor 
of exp(p 2 / 2 )- 

Turn next to the point (a) and (d) in the above list. From 
the biases in the relations /Aifw-^icns and ^i cns -^ n oisy found in the 
mock data, we derive an approximate mean relation of z^noisy — 
t^nfw = 0.07. We take this bias by simply shifting the counts as 
^noisy *" ^noisy (i/ + 0.07). 

Finally, the point (b) is corrected by simply increasing the 
peak counts by 10 percent. 

In Figure lBTl the corrected model prediction is plotted and is 
compared with the results from numerical experiments. To com- 
pute the corrected model prediction, n no i S y(V ), we first computed 
the uncorrected counts n pea k(i y ) at a small interval of 8u using 
eq. i\b\ . Then the local exponential slope of the counts is eval- 
uated by a finite differential scheme p{vi) = [logro(!/j+i) — 
log n(vij\/8v. We compute n no i ay (y) using the eq. IB2I with this 
p{vi). Then, the correction n no isy(f) — » n no i sy (/y + 0.07) was 
made for p(yi) at every Vi to correct the points (a) and (d). Finally 
the counts is increased by 10 percent to correct the point (b). The 
corrected prediction agrees with the noisy counts reasonably well 
at least v > 5, where the above approximations are valid. On lower 
peak heights v < 4, the prediction is systematically smaller than 
the measurement, where false peaks due to the noise have larger 
contribution. Note that the amplitude of the effects (b) and (d) may 
depend on the noise properties and on the choice of the smoothing 
scheme. 

It is clearly seen in the Figure IB 1 1 that the corrected peak 
counts has a slope very similar to the un-corrected counts. There- 
fore, one may make a simple correction to the model prediction ob- 
tained from eq. <16t by multiplying the factor 1.2, which gives an 
alternative approximate way to compute the corrected peak counts 
with a reasonable accuracy. 



APPENDIX C: RELATIONS BETWEEN THE HALO 
SHAPE AND ft AMPLITUDE 

Here we examine the effect of the individuality in the halo mass dis- 
tribution on the amplitude of the k peak. Figures lcT1andlC2lshow 
how the ratio lAens/fnfw varies with the halo shape parameters Q 
and R (see equations <20t and 1231 for the definitions), where we 
consider the class-(i) halos with f n f w > 4 and ^i cns > 4 only. The 
filled circles show the average over halos sitting in the range de- 
noted by the x-direction bar, where each interval is defined so that 
it contains 35 halos, and the error bar in y-direction denotes the 
RMS among them. 

The correlations between the shape parameters and amplitude 
of the lensing signal are observed in both the figures, though the 
statistical significance is not very high. The anti-correlation in the 
!4ens — Q is because more concentrated mass distribution (Q < 1) 
compared to the average NFW profile provides greater lensing sig- 
nal. The lensing signal is also amplified if mass distribution within 
a halo is elongated along the line-of-sight (R > 1). However, we 



© 2003 RAS, MNRAS 000.[Tll22l 



22 T. Hamana, M. Takada & N. Yoshida 




Figure CI. Shown is the scatter plot of the distribution of the halos in 
Q/QNFW-^lcns/'-'nfw plane. The class-(i) halos having z/ n f w > 4 and 
^lcns > 4 only are plotted. The filled circles show the average over 35 
halos within bins whose ranges are denoted by the x-error bars, while the 
y-error bars denote RMS among the 35 halos. 
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Figure C2. Same as Figure lcTl but for R—V\ en s / v niv/ plane. 



draw attention to the rather large scatters in these relations that may 
arise from the following two sources. First, the shape parameters 
we employ characterize rather global mass distribution within a 
halo and cannot describe the detailed structures. Second, probably 
more important one is the projection effect due to different struc- 
tures along the same line-of-sight of the halo. 
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